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ABSTRACT 


The numerical frontogenesis model of Williams (1973) 
is modified to include moisture with its subsequent con- 
densation and release of latent heating. The turbulent 
diffusions of momentum, heat and moisture are repre- 
sented with various coefficients. The numerical solutions 
show realistic quasi-steady fronts forming within one to 
two days. These solutions are Ae and compared over 
a range of the various coefficients, and various ranges of 
temperature. 

Inclusion of moisture in the model causes intensifi- 
cation of baroclinicity at mid and upper levels. AI1so 
noted is a sensitivity of the moist model to reference 
potential temperature due to the exponential relationship 


between saturation vapor pressure and temperature. 
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I. INTRODUCTION 


Hoskins and Bretherton (1972) have shown analytically, 
and Williams (1967, 1972) has shown numerically that dis- 
continuous fronts can form within a finite period of time 
if no turbulent ‘diffusion is present. These studies sug- 
gest that a discontinuity will form within 24 to 36 hours 
when reasonable initial conditions are used. If turbulent 
diffusion is present, it can be expected that a balance 
will be achieved between the frontogenetic advections and 
turbulent diffusions of momentum, heat, and moisture. The 
front should remain in this state of quasi-balance as long 
as the large scale deformation field causes frontogenetic 
advections around the front. 

The numerical frontogenesSis model of Williams (1973) 
shows the effects of horizontal and vertical turbulent 
diffusions of momentum and heat in the formation of these 
quasi-steady fronts. This model, which also includes an 
Ekman boundary layer, will hereafter be referred to as 
W73 and the model as it existed in Williams (1972) will be 
referred to as W72. 

This investigation appends moisture with its subse- 
quent condensation and release of latent heating to W73. 
It is assumed that all the latent heat released in this 


condensation is used to heat the air isobarically. 





The purpose of this study is to obtain and examine 
frontal solutions from the model throughout a four day 
period. The physical model is essentially the same as 
one of the models treated by Hoskins (1971) and Hoskins 
and Bretherton (1972). The basic model, as it now exists, 
numerically describes frontogenesis which is forced by a 
nondivergent horizontal wind field which contains stretch- 
ing deformation. This deformation wind field is constant 
in time and is independent of height, except in the Ekman 
boundary layer. The moist hydrostatic primitive equations 
with diffusion are used in which the time dependent quan- 
tities are functions of Y and P only. Therefore, output 
from the model is in the form of vertical cross-sections 
in the Y-P plane. | 

W73 has been transformed to pressure coordinates to 
obtain more accurate values of pressure for the moisture 
equations that have been included in the present model. 

In Section 2, the basic forecast equations are devel- 
oped, transformed to pressure eooLaeees and a simplifi- 
cation which keeps the problem two dimensional is discussed, 
The initial conditions are presented and transformed to 
pressure coordinates in Section 3. Convective adjustment 
processes and latent heat release in the model are included 
in Section 4. The experiments undertaken are examined in 


Section 5 and conclusions are given in Section 6. 





II. BASIC EQUATIONS AND BOUNDARY CONDITIONS 


The hydrostatic primitive equations with diffusion 


and moisture May be written as follows: 











av d (wv) ; a°V 
eee orm An! SPSa TD eo 
a8 og. 3 7 2 976 
— VNC OW) ze (WO) = Ag O+C, m2 Teo), + Chae, C22) 
oq +V°(qv)+ oe ai A V7q+c g + M -—- M C2ie,) 
ot ~ ds q qd 5.2 A c : 

= - ds E : 
90> = L(s)98 ; (2.5) 


where the vertical coordinate has been replaced by s to 
aid in later discussions on transformation. The equations 
are from W73 with slight modification due to the addition 
of the equation of conservation of water wapor (2.3) which 
is analogous to the thermodynamic equation (2.2). The 
function L(s) in the hydrostatic equation (2.5) is a con- 
stant g/8 | in height coordinates. Also sinks and sources 
of heating and moisture have been added to the thermody- 
namic equation (2.2) and the conservation of water vapor 
equation (2.3). 

The notation used is basically the same as W72 which 


uses height as the vertical coordinate. However, in our 
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subsequent transformation to pressure coordinates, § is 
defined as a departure from the domain average of the 
initial state, $¢? represents geopotential, gz and w is 
ds/fdt in the continuity equation (2.4), thereby allowing 
compressibility. 

The quantities, An! Ags and at: are the horizontal ture- 
bulent diffusion coefficients of momentum, heat, and 
moisture, respectively, while Cat Cor and aa are the cor- 
responding vertical coefficients. These constant coeffi- 
cients were assigned various reasonable values in the 
different metist une experiments. Dry and moist convective 
adjustment processes are represented by the functions Qn 
and Qu while latent heating due to large scale condensa- 
tion is represented by Hae In the moisture equation, My 
is a Sink or source of moisture due to moist convective 
adjustment while Ne is a sink for moisture due to large 
scale condensation. 

Boundary conditions for the above equations are such 
that the domain is bounded by two rigid quasi-horizontal 
Planes. The no-slip condition at the lower boundary and 
the no-stress condition at the top roughly represents the 


troposphere. The conditions are as follows: 


WirevenO,t) = w(x,y,H,t) = O 
V(x,ry,oO,t) = 0 
ay | 
oe Ze (Xe YeH,t) = 0 7 Zo) 
089 e 06 a 
Cy ze l*r¥,0,t) = Co ye (%rV eH, t) = 0 
og _ dq 
5 Ze (XV Ort) = ee Jeg (Xe ¥eH,t) = QO 


Jane 





eee H represents the distance between the upper and 
lower boundaries. The lower boundary conditions are realis- 
tic for fronts over level terrain. Therefore, fronto- 
genesis in the lower troposphere is more realistic in the’ 
model than are conditions near the “tropopause", The addi- 
tional boundary conditions for moisture are completely 
analogous to the heating conditions in that fluxes of the 
quantities do not pass through the boundaries. 

The approximate steady state solutions to the hydro- 


Static primitive equations from [(2.1)-(2.5)] are: 


V=U= niles” Cos Bs) + onan) Bs] i 
-~ Iy(i=e “cos Bs)- ee etn Bs] jl} 
w= 0, ’ (2.7) 


¢?6 = 9 = amen + ane - fD xy + F(s) 


oe— G(s) 
q = H(s) 

eye - 
where B = (£/2C_) , D and £ are constant and where F(s) 


is such that the hydrostatic equation (2.5) is satisfied. 
The functions G(s) and H(s) are arbitrary functions of s 
Since w = 0. These solutions satisfy the boundary condi- 
tions (2.6) and allow for inclusion of the Ekman boundary 
layer. These relations (2.7) satisfy the steady-state equa- 
tions if some small advection terms in the boundary layer 
are neglected. These small terms are uSually dropped in 


the Ekman Theory. 


ee 





The equations and the model become two dimensional 


after a subdivision of the dependent variables as follows: 


V = U(x,y,x) + uly,S,t)i + Viy,s,t)j 

w = wl(y,s,t) 

Boe wey esc) : (2.8) 
¢ = @(x,y) + Tly,s,t) 


q = qly,s,t) 


It is noted that all departures from the approximate 
steady deformation solutions (2.7) are assumed to be in- 
dependent of x. If we substitute the relations (2.8) in 


the i component of (2.1) we obtain: 








du 3 0 -Bs 
— + =~ + = + - 
a+ Ny (uv) ye (wu) uD(l-e cos Bs) 
=Diviii-e "cos BS)oRe asan Bs 15 + vDe ?*sin Bs 
2 
+wD ae eo (eos Bs)+ye eT Bs]=fvt+A = 
os m 2 
dy 
2 
+c. qu (2.9) 
957 


Note that the assumption that u is independent of x 
for all time is violated in (2.9) in the two terms where 
X appears. This indicates that the other quantities 
which were assumed independent of x are, in fact, x- 
dependent and this would bring other terms into the equa- 
tions. In this study we will neglect these added terms 
and apply (2.9) and the other equations at x = 0. If the 


initial u,v,w,9,71 and q fields are independent of x, the 


3 





error from this approximation will grow slowly and will 
be confined to the boundary layer since the x-dependent 
terms in (2.9) are zero outside the boundary layer. 

Since frontogenesSis occurs very rapidly it is expected 
that the development of x-variations in the dependent 
variables would have only a small effect on the resulting 
quasi-steady front. In any case these effects could not 
be observed in the atmosphere, because atmospheric fronts 
have some variation in the basic fields along the front 


which would be much more important. 

















When equation (2.9) iS evaluated at x = 0 it becomes 
du d (uv) d (wu) oV du voV. -Bs . 
— pein oe ane ae Nees ce 
Beat ay + ys a ae ee sin Bs) 
2 2 
-~Vw ie Saeed Bs)=fvt+A gu + C me ; : (22 16) 
ds m 2 m 2 
Jy ds 
where V = -Dy and [. = eeesecs Bs.e- The j} component of 
(2.1) applied at x = 0 takes the form 
dV 0 d 0 oV -BS., 
re + as eer a og e Sin Bs 
2 2 
reer Cosmas =<" -cume® 2 +e 2. (eit) 
ds dy m dy? M 952 


In a similar manner from equation (2.2), the following 
equation for the departure of potential temperature can be 


derived: 








2 2: 
3° i ates) r —_ ry ge z A 0“6 ip ee a°6 
y) S y dy? 52 

a oe cf en + He ’ (2722) 
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and the following equation for q follows from (2.3), 








2 2 

oo, 20a) . S-wwq) + tv 2 -a 294 ¢ Fo 

y oy g dy? g ds2 

Equations (2.4) and (2.5) become 

OV Ow | 
ee Ss (24an 
and 
om 
ee Soe (2.15) 


The boundary conditions at the bottom and top of the domain 


respectively are 


ae a0 _ oq _ 7 
du _ dV _ _ 00 _ dq _ = 
ee on vee Cg os ae as ot = (2.16) 


If we define the vertical average of a quantity as 


H 
f (¢( )ds, (2217) 
O 


A 
Mv 
tl 

mir 


integrate the hydrostatic equation (2.15) with respect to 
Ss, and remove the vertical mean, it follows that 
. Ss S 
TM-<tr> = [ J L(s)8ds- < f L(s)9ds>] . (2.18) 
oO O 


Take the vertical average of (2.14) and use the boundary 


(2.16) which gives 


— <v> = 0 , (2.19) 
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This equation states that the total mass flux in the 3 
direction is independent of ys. In W72 a symmetry argument 
was used to show that this flux must be zero. This argu- 
ment does not hold strictly in the present experiment, 

but it can be expected that at large distances from the 
frontal zone the disturbance mass flux will vanish. Thus, 


it is assumed that 
3 (24.2.0) 


If the development of W72 is used, equation (251) can 


be written 


oy 3 0 3 
—_——= Eee = —— Se = 
rE + pe <vv>Jt+ ye (WV) + aye <Tv> V) 


OV -Bs . -Bs_. 3 -BsS 
- ~[ue sin Bs-<ue sin Bs> -Vi{w ~-(e cos Bs) 


oy os 





- <w a BS)>] = - OG =<m>) -f£({u-<u>) 
os dy 
2 2 
7 ae c 1 aan =(2%) 4) “ ; (2.21) 
dy? ds7 | 


The equations (2.10), (2512) "C2.13), (2. 84), (27s) 
and (2.21) form a complete set which can be solved by a 
pure marching process. The finite difference equations 
conserve mean squares in the advection terms and are de- 
scribed by Williams (1967). The Matsuno scheme is used 
everywhere except in the arrangement of variables and the 
finite difference approximations are the same as those 


used by Williams (1967). In order to close the problem, 


16 





computational boundaries must be introduced in y. Since 
the disturbance velocities should die out at a sufficient 


distance from the axis of dilation, then 


v(+ Y,s,t) = 0. (e522) 


However, there is appreciable inflow across these com- 
putational boundaries since V(+ ¥) = + DY . The quantities 
u, 98, and q which are advected across the boundaries must 
be specified independent of the interior values if compu- 
tational stability is to be maintained (Platzman, 1954). 


Thus, the following boundary conditions are used: 


Usiet CY oc Bye 2), Ss; & u PCy, + By72) 7S, 0)! 


O[+(Y + Lge?) St) SOY ny 72), 5,01 ° (2.23) 


q[+(¥ + Ay/2),s,t] q [+(¥ + Ay/2),s,0] 


The computational boundaries y = + Y are placed between 
gridpoints so that the above conditions are actually 
Sppliedmeat y jet. (Y + fy/2). 

The boundary conditions were found to be satisfactory 
except that a solution separation developed near the 
boundaries. This was controlled by using a Matsuno time 
differencing scheme for stability. 

In order to have precise values of pressure at grid- 
points for the moisture experiments it was deemed neces- 
sary to transform W72 to pressure coordinates. 


Minimization of changes to finite differencing schemes 


i 7 





of the model have been achieved by defining a parameter gs 


which is analogous to a field of z. Specifically, 
s = P ar og i (2.24) 


a field that will increase in the upward direction and 
from which values of pressure can easily be obtained. 
Further, it is seen that ds = - OP so that finite differ- 
encing schemes in W72 can be equally well applied to the 
field s. 

Through the expression in Ss coordinates the momentum 
equations (2.10) and (2.21), the thermodynamic equation 
(2.12) and the moisture equation (2.13) remain as before 
for the model. However, the continuity equation (2.14) 
in the model actually remains in w and describes a nega- 


tive field of omega in that, 


= 
It 
eal 
ri]M 
It 


ee (2525) 


Equation (2.5), the hydrostatic equation, differs 
Significantly due to the transformation. The derivation 
is briefly as follows. The hydrostatic equation expressed 


in pressure coordinates and the equation of state gives 


ee eg es Sig (Zi) 


If a form of the potential temperature definition is 


used, 
r= o(Py)* , eT) 
Po 7 
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and combined with (2.26) and (2.27) it follows 


9d _ _ 
dp 


oO |w 


(B-)" 9@., (2528) 
Po 


It can be shown that 


R K 3 K 
= (=) = ¢ =) (2.29) 
PT POP Ree 
The latter form simplifies finite difference energy conser- 
vation. Therefore, substituting equation (2.29) into 


equation (2.28) and transforming to s coordinates we obtain 


se Pee 1h) ") 6 , (2.30) 
S 


a form of the hydrostatic equation. Note that L(s) is not 
constant as in height coordinates since 


ee cence (2.31) 


Po p 9s Pp. 


L(s) = 


'O lw 


It follows that the complete set of equations solved by 
the marching process remains intact but with L(s) defined 
acino( 2.31). 

The boundary conditions (2.16) also remain intact 
except that the lower boundary of vertical motion in the 


model is approximated as 


w=-wWwWe= 0,ates=oaQ., (2232) 


io 





Iii. INITIAL CONDITIONS 


The initial conditions in height of s coordinates are 
the same as those employed in W73. The initial potential 


temperature field is given by 





96, H 2 
O8(y,z,0) = (— (z- 3) ce (=) arctan (Sinh @ y) 
where 
Q = ao a i (3.1) 
06 
ir 
HWL(s) 2 


The quantity 96, /9z Wwimilen ene COnSstant is the initial statue 
stability, or being the initial horizontally averaged 
initial theta and "a" is one-half the total horizontal tem- 
perature variation. 

By substituting equation (3.1) into the thermal wind 
equation and integrating we obtain the x~component of the 


rotational wind which is 


a (z - >) sech GQ y= (3.2) 
oO 





WG Y¥. p20) = 


A] 


The divergent part of the wind or v field is obtained 
from the quasi~geostrophic equations in z coordinates. 
Details are as in Williams (1972). The remaining fields 
of w and 7 are obtained from the continuity (2.14) and 


hydrostatic (2.15) equations. 
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Once again transformation of the above initial fields 
must be made to pressure coordinates. 

To obtain initial conditions for the theta field we 
must first find 8, , the horizontally averaged initial 
theta based on O, (po) - Derivation is as follows. Given 


knowledge of differentials and the definition of geopoten- 


tial it is seen that 





mt oe a Or 2g (30a) 
dp dz op g 0z dp ’ i 


where 90 /dz is the prescribed static ‘stability at t = 0. 
Substitute a definition of potential temperature (2.27) 


to obtain 


90 30 a 
= eee ele (3.4) 


0 OZ 





Wilke 
to |x 
ay 


6) 90 
fe spear f - Esthet ae, cae 
6 p g 2 Ee Po 


= (Se (3.6) 
Pp 





where K = a , it follows that 


7 < Cc. 930 


es iP K 
O, (Pp) 0, (po) EXP [o a ) 





yyy. (3573) 
p 


O 
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oe =i wh 








: 
- om o————=» =& & @ 


—————— oe ee 
= 
i i) 


“4 


Equation (3.7) is representative of term one of equa- 
tion (3.1). Therefore, the initial potential temperature 


field in pressure coordinates is given by 





5 cp 21, eye 
O.(y»p,o) = 0, (p.) Ble Te on! )] 
-a(2/m) arctan (Sinh @ y) . (3.8) 
where 
a 2 el al (3.9) 









30, 
L(s_) ae 





In the model @ is given a constant median value by 
solving equation (3.10) at p = Dea 60 cb and using the 
static stability 90,/9p at this level. Recall also that 
pressure in equation (3.8) can easily be obtained for 
model purposes as p = Po 7 S- 

Further it should be recalled that @ is a departure 
from the domain average of OF. This domain average is 
found in the model by finding the vertical average of 
equation (3.7), < Cee 

The initial x-component of the velocity or initial 


u equation is found using a form of the hydrostatic equa- 


PlOnm2.2o)eana the definition of geostrophic wind, 
a 
UC he a Tee ? (3.10) 


ou _ 
dp 


—_ = 


bh] 


Pee Seep kK 00 
ae Pas . (ao0t) 


ue) 
QoQ 
het 
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If we again assume a median pressure level Dees 60 -cb, 


and that at this level u(p) = 0, we integrate 
Pp Pp 
/ = dp -< = or dp . (areiey 
Pn Pn 
Recalling that u(p) = 0 and that K = R/C, the integra- 


tion yields 


C 
u(y,p,0) = => Saag P sie (= a) (3.13) 


dy Po re 


From (Williams, 1972) equation (AQ3) ae asl seenmenat 


06 _  2aa 





sech (ay) . (3.14) 


Therefore, the initial u field in pressure coordinates is 


given by 


” 


p 
is Faia wi em) sech (ay) . (euisp 
© © 


2aad 


u(y,p,o) =e 7 


The initial divergent wind v is determined as in Wil- 
liams (1972) but using the quasi-geostrophic equation in 
pressure coordinates. The stream function is defined as 

v = d¥/dp 


’ (35 lo) 
and Ww =-dp/dy 


and equation (4.10) of Williams (1972) takes the form 








2 36 2 
5 a L(s_) (- = ee £2 Lis} o (3.17) 
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where L(s) is assumed 


The equation (3.17) is solved for Ww from the initial 
temperature field (3.8) with the technique of Ogura and 
Charney (1962) and v is then determined from (3.16). 

It is noted that these initial conditions do not 
Satisfy the boundary conditions on u, v, and 8 at the 
upper and lower boundaries. A period of adjustment will 
be required to form the eer ereeion layer and some 
oscillations may be observed later. 

The fields of w and 7 are obtained through the con- 
tinuity and hydrostatic equations (2.14) and (2.15). 

Finally it is necessary to define explicitly methods 
used in determining the initial field of specific humidity, 
given O, (y,P,0) from equation (3.8). 

Saturation vapor pressure Cur is first determined at 
each gridpoint through use of an integrated form of the 
Clausius~Clapeyron equation, an equation expressing satura- 
tion vapor pressure as a function of temperature. Accord-~ 
ing to Hess (1959) 


Tye: 


. ae 
e. = .611 EXP[ e ae (3.18) 


a 
R 
v 


where e. is expressed in centibars, Ly is the latent heat 
of condensation, RY is the gas constant for water vapor, 


and TO is a constant of integration. 
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Knowing saturation vapor pressure, the saturated speci- 
fic humidity was obtained from the approximated relation 


from Haltiner and Martin (1957) 


-622e 
S 


a mon ° (Sie 


Therefore, by having knowledge of potential temperature 
at t = 0 from equation (3.8), q, can be found at each point 
uSing equations (3.18) and (3.19). By defining q to be a 
given fraction of qe in different experiments, atmospheres 
of different constant initial relative humidities were 


tested. 
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IV. CONVECTIVE ADJUSTMENT, LATENT 
HEATING AND PRECIPITATION 

Inclusion of two additional heating terms in the moist 
thermodynamic equation (2.3) necessitated some further 
amendment to W73. Atmospheric heating by dry convective 
adjustment, Qn: was included in W73 to eliminate any in- 
tense grid scale convection which developed due to the for- 
Mation of any dry unstable lapse rates. This redistribution 
of heating was accomplished uSing averaging techniques to 
ensure that the dry adiabatic lapse rate was not exceeded 
3. < 0). 

Lapse rates were tested in each vertical column of the 


09 


model. Upon encountering an unstable layer aa > O)}>3 a 


cumulative average was taken and applied to the gridpoints 
within this layer until neutral stability (se = 0) was 
attained. Simultaneously the model was checked for insta- 
bilities developing below the convectively adjusted layer 
and potential temperatures at these levels were also cumu- 
latively averaged and reassigned the new values. The final 
result of the technique is that the heating is redistri- 
buted to simulate physical processes of the atmosphere and 
mathematically that + SUS ehrougnhout cach colunn, of, the 


model. A completely analogous scheme was used for saturated 


parcels in the moist experiment, only uSing equivalent 
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potential temperatures, =, in place of potential tem- 
perature, 6. 

With the addition of moisture,consideration of large 
scale condensation and the subsequent release of latent 
heating became important. 

According to Haltiner ( ), changes in moisture and 
temperature due to large scale condensation were deter- 
mined as follows. Where relative humidity exceeds 100%, 
the adjustments in temperature 6T and specific humidity 


6q are determined by the following equations: 


q + dq q(T + 6T, p) , (Aes) 
and 


ClOTe=- Legoaqe. 4.2 
P q ( ) 


These equations require that the final state of the 
air be exactly saturated and that excess moisture be con- 
densed isobarically, releasing the latent heat to the air. 

Equations (4.1) and (4.2) can be eee in itera- 


tive form as 


re 


Ss 
aT ys AT ? g a Ee ? (4.3) 





Aq = q, (1l-r) + af 


and Cel = ie ENG) . (4.4) 


Eliminating AT from (4.3) and (4.4) leads to 
qe (Ly) 
Aq SS = ° 2) 


dq 
L Ss 
tte GP? 





od 


=@ 


i = iba = SB 


—.e ee 





At constant pressure the Clausius Clapeyron equation 


can be expressed as 





dq iad 
(a = ——, (4.6) 
Pp Re 7 


Substituting (4.6) into (4.5) we obtain 


q_(1-r) 
Aq == .. (4.7) 


Now Aq may be calculated from (4.7) and the correspond- 
ing value of AT from (4.4). 
The first estimates of the adjusted temperature and 


Specific humidity are, therefore, Se 
ie eter and q° = sq + Aq . (4.8) 


There values in turn were used to obtain improved ap- 
proximations. Iterations were containued until the values 
of T" and q‘' at each level matched with sufficient accu- 
racy one hundred percent relative humidity. The changes 
were then applied to the potential temperature and speci- 
fic humidity fields of the model accounting for the term 
H A Of equation (2.2) and M of equation (2.3) 

The model was then investigated for unstable moist 
convection and its associated heat of condensation. A 


convective adjustment scheme similar to the dry one was 


used on the saturated parcels so that moist static stability 
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was always maintained. A field of equivalent potential 


temperature was formulated according to Holton (1972) as 





(4.9) 


and the moist convective adjustment was designed to ensure 
thateehbesmorst saddabatic lapse rate was not éxceeded 
(00 _,/9p < 0). 

The heat release, Qn? as a result of moist convective 
adjustment was included through knowledge of the change of 
equivalent potential temperature during the adjustment. 

By taking the natural logarithm of equation (4.9) it 


is found that 


ln an = In O + eer 





(4.10) 


In a model using pressure coordinates, a constant pres- 
sure surface may be assumed, so by differentiating both 
sides of equation (4,10) and recalling the definition of 


potential temperature gives 


dO» 7 aT a oe ae = ne qs OT (4.11) 
2 a e 
oS. sy c. 7 “ 


Solve the Clausius Clapeyron equation (4.6) for dqe. 


Substitute in equation (4.11), and solve for dT. The 
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solution is found to be 











do: 
at = 2 
q_L Lq 
1 s Cc S 
Oe [FQ + > C =! (4.12) 
RCT p 
Vp 
. 90» 
where d0_» is the change in aan needed to adjust ap = 0. 


Since the atmosphere remains saturated throughout this 
process of moist convective adjustment, the Clausius- 
Clapeyron equation was again used to obtain the ee. sea 
ing dq or actually dq Once again the changes were then 
applied to the potential temperature and specific humidity 
fields of the model accounting for the term Qn of equation 
(2.2) and M Of equation (2.3). 

A cumulative amount of precipitation was then found 
using the relation 

Po . 
p=- : 6q dp (4.13) 


oO 
in agreement with the changes in specific humidity due to 


condensation and moist convective adjustment. Precipita- 


tion was expressed in inches of rainfall in the model. 
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Ve NUMERICAL SOLUTIONS 


All numerical results to be shown use the following 


values for the constant: 





At = 540 sec Po = 100 cb ene 60 cb H = 80 cb 
Ap = 1.78 cb Y = 1800 km Ay = 60 km g = 9.81 Mm sec. 
f = syle D= TO" S60" se 300° 
O,(p,) = 270, 280, 290, 300°K C. = WOGEM Ou lest aoe 
R, = 461 joules oe Oe R = 287 joules es 
6 = = 
K = R/C L = 2.5 x 10 joules kg . TT. = 273°K 
p Cc Oo 
a =12.56°K RH = 50,70,85% 
a0 30 
a 4°K km i ; (5) = -.53448, .6, .~8°K cb 7 

4 dp Po 
A =A, = A_ = 3 x tote 10> ii; ee 
m 8 q 
Ce = 0, . 0021409, - 30001761 a (which corresponds to 

ORetO ; bhsezs “eee in height coordinates) 
Cy = 0, .~001409, .0001761 priddeces (which corresponds to 

OO, S25 Meee height coordinates) 

2 = il ; 2 ae 

Ca = 0, cb sec (which corresponds to Om sec in height 


coordinates) 
where more than one indicated value means eiait the constant 
was varied in different experiments. Whenever possible, 
the values are the same as those used in W73. 

In this section numerical frontogenesis solutions with 
moisture are examined at t = 4 days and compared to the 


corresponding dry atmosphere. 
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As can be seen from Table I, 26 experiments were per- 
formed on five cases, cases selected from Williams (1973). 
Case 1 contains only horizontal turbulent diffusions of 
momentum, heat, and moisture while cases 2 through 5 con- 
tain both horizontal and vertical turbulent diffusions 
of these parameters. Values chosen for these turbulent 
diffusion coefficients were computed to correspond with 
values used in Williams (1973). 

In Table I, 8, (po) is the reference potential tempera- 
ture (°K) as introduced in equation (3.7). The condensa- 
tion interval indicates the number of time steps between 
Successive applications of condensation and moist convec- 
ee adjustment to the model. Turbulent diffusion coeffi- 
cients are then given as previously defined followed by 
the initial relative humidity (percent) and the cagieiad 
Static stability (°K/cb). 

Experiments were first performed with a reference 
potential temperature (8, (po)) or S00 Keto ailrectly compare 
with the dry solutions of W73. This corresponds to a tem- 
perature range of 62-104°F at the surface along a North- 
South axis. Therefore, experiments were also conducted 
with Op» = 270,280,290°K to simulate physically realis- 
tic situations. Changes to the dry atmosphere with the 
different references were almost negligible, but due to 
the exponential function in the Clausius-Clapeyron equa- 
tion important differences were observed in the moist 


atmosphere. 
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Experiments were also conducted using different values 
for initial relative humidities and initial static stabili- 
ties. 

A brief tabulation of experimental results is given in 
Table II, comparing widths of frontal zone, frontal slopes, 
vertical motions and precipitation. 

A reasonable: measure of the width of the frontal zone 


is given by 


ac O0(-Y,s,0o) - 6(Y,s,0) 
a 1.28 f 


Oy ‘max 


Goi) 


where so is approximated by a one-sided difference. This 
width (km) is indicated at three levels for each experiment 


in the table, 
a6 


ae ee at three levels is 


The North-South location of | 
y, (km) and is used aS a position indicator of the frontal 
zone where ly,| < 1800 km. 

The pressure level and magnitude of maximum vertical 
motion is then given followed by the North-South position 
of maximum precipitation (km), ly, < 1800 km, and the 
corresponding amounts of cumulative precipitation over the 
four day period. 

It is noted from Table II that widths of frontal zones 
between corresponding moist and dry experiments are pri- 
marily effected at mid to upper levels with only small 


deviations being noted at the surface. The expected in- 


crease of thermal gradient at these mid to upper levels 
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due to latent heating in the moist cases is seen to be 
the cause. 

Comparisons of frontal slopes of corresponding experi- 
ments revealed deviations again at mid levels. The slopes 
appear more vertical throughout these levels in the moist 
experiments. Graphical illustrations will later be pre- 
sented. 

Maxima in vertical motion saen to intensify and be 
found at higher levels along the frontal zone in the cor- 
responding moist cases. This effect can also be directly 
tied to releases of latent heating in this region. 

commas precipitation amounts in Table II are ob- 
served to be primarily functions of initial relative 
humidity and reference potential temperature. 

Primary emphasis on analyses of experimental results 
waS on contour plots from the IBM 360/67 computer. These 
results obtained for all 26 experiments are analyses at 
t = 4 days of the fields of u,q,8' and w where @6' = 6 - oF 
and 6, is the horizontally averaged initial potential tem- 
perature. The entire vertical plane is shown from 100 to 
20 cb in these plots. However, the outer portion of the 
domain between ly | = 1200 km and ly | = 1800 km is not shown. 
Also note that isolines of specific humidity in the figures 
to follow have different intervals dependent upon reference 


potential temperatures in order to allow for proper defini- 


teonw@ofethe field. 
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Due to the large number of contour plots obtained, ex- 
periments 1,2,4,5,6,11,12,13,17,18 were chosen as repre- 
sentative samples for detailed discussion. Experiments 
1,2,4,5,6 and 11 contain only horizontal turbulent diffu- 
sions of momentum, heat, and moisture while experiments 
12,13,17 and 18 contain both horizontal and vertical turbu- 
lent diffusions ‘of these parameters. 

Figure l comparing Sporimenus l1 and 2 shows the 
frontal scale at 99.1 cb reaching steady-state at approxi- 
mately t = 2 days. Little deviation is seen between the 
wet and dry experiment. Figure 2 compares the frontal 
slopes of the same two experiments. AS mentioned before 
Slopes are seen to be largely the same with the only signi- 
ficant deviation being at midlevels where the moist slopes 
approaches the vertical and then recurves in agreement 
with the dry solution. 

Again this deviation in the wet experiment is probably 
due to the release of latent heating at these levels. 
Similar results were achieved for other cases with more 
extreme deviations being observed at higher reference 
potential temperatures. 

Contour plots for experiments 1 and 2 are shown as 
figures 5 - 12. Figures 5 and 6 show distributions of the 
rotational wind, u for a corresponding dry and moist experi- 
ment respectively. It is noted that low level features are 


guite similar in the easterly flow. However, horizontal 
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and vertical wind shear iS more intense at low mid to mid 
levels in the warm air in the moist experiment. It is note- 
worthy that this property is dominant throughout the moist 
solution. 

Figures 7 and 8 contrast values of @' throughout the 
domain. Noted is the more intense gradient of potential 
temperature through the mid levels and the pocket of warm 
air at mid levels most likely due to release of latent 
heating in the moist experiment. 

Fields of specific humidity are contrasted for experi- 
ments 1 and 2 in figures 9 and 10. In the dry experiment 
moisture is merely advected and horizontally diffused 
throughout the domain showing an interesting though physi- 
cally unrealistic situation. The processes of condensation 
and moist convective adjustment are only applied in the 
moist experiments. By comparing the isolines of specific 
humidity, it is seen that the processes of condensation 
and convective adjustment are indeed reducing the moisture 
content of the air and that the gradient of specific 
humidity is stronger along the frontal zone in the dry 
experiment where no condensation is taking place. 

Vertical motions are investigated for the experiments 
in figures 11 and 12. Upward vertical motions are stronger 
and occur at higher levels closer to the frontal zone in 
the warm air mass for the moist experiment. Downward ver- 
tical motion in the cold air mass was also slightly inten- 


sified with its maximum also at a higher level. 
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Figures 13 - 28 compare the dry and moist fields of 
u, gq, 9° and w at reference potential temperatures of 
290, 300 and 270°K. Figures 13 ~ 16 show moist fields 
(8, (po) = 290°K). The horizontal and vertical shear of 
uin figure 13 is seen to be stronger and to be extending 
to higher levels than in the previous moist experiment. 
Stronger thermal gradients are-Seen to extend higher in 
the atmosphere in figure 14. Latent heat release occurs 
at higher levels and to greater extent as is seen in the 
upper level isentropic pattern to the warm air side of 
the frontal zone. 

IG is predicted by the ClauSius-Clapeyron equation 
(3.18), the moisture values (figure 15) are significantly 
greater than in the previous experiments. Again note 
the decreased gradient of moisture along the frontal 
zone in the moist experiment. 

Vertical motions, shown in figure 16, are also somewhat 
intensified, occur at higher levels and have slightly 
more intense gradients along the discontinuity in the warm 
wie. The downward vertical motions in the cold air are 
also somewhat increased, 

Figures 17 ~- 24 contrast dry and moist experiments 
with a reference potential temperature of 300°K. Little 
difference is seen between the dry experiment and the 
earlier dry fields based on 280°K except for a northward 
displacement of the fields at the colder reference poten- 


tial temperature. Due to the exponential nature of the 
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moisture parameter, moist fields showed significant changes 
in the following experiments. 

The rotational wind u, as is shown in figures 17 and 
18, has much stronger horizontal and vertical shear in the 
warm air at low mid to upper levels than in previous moist 
‘experiments and even has an area of easterly winds at upper 
levels in the warm air. The air 1S much warmer in the 
upper levels due to latent heat release, as indicated by 
figures 19 and 20. Thermal gradients along the frontal 
zone in the warm air are also stronger in the mid and upper 
levels. | 

Gradients in the moisture field are greatly intensified 
and specific humidities are in excess of 28 g/kg at low 
levels to the South. Recall that the "dry" q field (figure 
21) is distributed due to advective and horizontal diffusive 
processes only. Perturbations in the moist experiment 
(figure 22) are probably due to the larger quantities of 
latent heat release observed in this experiment. 

Vertical motions (figures 23 and 24) exceed 4 x oie 
mb/sec in the warm moist air. These features occur at 
upper levels and a cell containing downward motion is 
observed at these upper levels to the South and seems to be 
somewhat due to the impermeable upper boundary. 

Fields with a reference potential temperature of 270°K 
are finally examined in figures 25 - 28. Due to minimal 


moisture content at these lower temperatures the fields 


are quite comparable to dry experiments previously shown. 
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Shear in the u field (figure 25) and thermal character- 
istics (figure 26) are also quite similar to dry experi- 
ments. Specific humidity (figure 27) reaches a maximum 
of 5 g/kg allowing for only a small possible latent heat 
release. Pinally, vertical motions (figure 28) occur 
Mainly at the lower levels and are mcre closely tied to 
dry experiments. 

Further experiments with case l are included in 
Table I. Due to the repetitive nature of the contour plots 
obtained, they are not included as figures, but some fur- 
ther discussion is appropriate at this time. 

Experiments 3 and 8 were run with varying initial rela- 
tive humidities. Higher initial humidities do lead to 
more latent heat Bene with its increased thermal gradi- 
ents, wind shear, and vertical motion in the warm air but 
did not seem to warrant further experimentation. 

Experiments 9 and 10 were run with varying static sta- 
bilities but due to extensive similarities in the solu- 
tions, further experiments were not required. 

The condensation interval was generally five time steps. 
The iterative procedure of returning relative humidities 
to 100% worked equally well at every ten time steps, never 
exceeding two iterations. However, five time steps were 
finally chosen to compromise computer time and physical 
reality. Only minor changes were observed in the contour 


pEobSMoteusrq,0" and w. 
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Experiments that follow contain both horizontal and 
vertical turbulent diffusions of momentum, heat, and 
moisture. As in the horizontal case previously discussed, 
values of the diffusion coefficients were computed to 
correspond with those used in Williams (1973). 

Figure 3 compares the frontal slopes for the corres- 
ponding dry and moist experiments 12 and 13. Once again, 
Slopes are almost identical with the exception of at mid- 
levels, where in the moist experiment maximum potential 
temperature gradient again shifts toward the warm air and 
then recurves to agree with the dry experiment. 

Figure 4 for the same experiments shows frontal scale 
at 99.1 cb versus time with steady-state being obtained 
at approximately.t = 2 days. Oscillations arise from the 
development of the surface boundary layer, which is not 
present in the initial conditions as was described in 
Williams (1973). Very little difference is observed between 
the moist and dry solutions at low levels aren reinforces 
the idea that the model remains essentially unchanged at 
these levels with addition of moisture, 

Contour plots are then shown as figures 29 - 36 for 
experiments 12 and 13. Compared to the corresponding dry 
case, horizontal and vertical wind shear is again increased 
especially at low mid to upper mid levels in the warm air 
(figures 29 and 30). The thermal pattern (figures 31 and 


32) shows increased gradient at the low mid to upper levels 
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in the warm air and the area of maximum latent heat release 
is seen in the isentropic pattern. Vertical influences 

for the moist experiment seem Somewhat More pronounced for 

a reference potential temperature of 280°K as compared with 
similar experiments containing only horizontal diffusion. 

Figures 33 and 34 compare the fields of specific 
humidity. The changes continue to be a reduction in 
gradient along the frontal zone in the prime area of latent 
heat release in the moist eee Little change is 
observed throughout the rest of the domain. 

Vertical motion fields (figures 35 and 36) for these 
Same experiments show low level maxima of 1.5 x lige ey see 
in the warm air but an additional more intense maximum of 
vertical motion is found at mid and mid upper levels along 
the frontal zone in the moist experiment. 

Case 3 includes a vertical thermal turbulent diffusion 
coefficient eight times as large as in case 2 previously 
discussed. From Williams (1973) a decrease in-:frontal tilt 
is expected along with an increase of frontal scale due to 
the change in diffusion coefficients. This can be observed 
by comparing appropriate figures of the dry atmosphere. 
Other changes to the fields will most likely be due to in- 
clusion of moisture. At any rate our comparisons will 
remain as before, a comparison between the corresponding 
dry and moist experiments. 

Figures 37 and 38, comparing experiments 17 and 18, 


show again the increase of horizontal and vertical wind 
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Shear at low mid to upper mid levels to slightly greater 
extent than in the previous moist experiment, experiment 
30. 

The isentropic analyses (figures 39 and 40) show the 
tightening of potential temperature gradient at low mid 
to upper levels once again. However, it is now observed 
that the increased thermal diffusion has smoothed somewhat 
the perturbation in the isentropes in the area of maximum 
latent heat release. This added vertical thermal diffusion 
seems to aid in distributing the isolated heating induced 
“pocket" that had been noted previously. 

Figures 41 and 42 compare specific humidity fields of 
the same experiments showing the weakening of gradients at 
low mid to upper levels due to condensation. 

The fields of vertical motion, figures 43 and 44, are 
Similar in results to the previously discussed experiment 
containing horizontal and vertical turbulent diffusions 
except for some smoothing of features due to the added 
vertical thermal diffusion of the present experiment. Once 
again the maxima noted in the dry experiment 18 is main- 
tained with an added maximum in vertical motion in the mid 
level warm air along the frontal zone in the moist experi- 
ment. 

After an analysis of the WerOre mentioned experiments, 
it becomes obvious that an observational pattern has 
developed and that further discussion using the line of 


attack would be largely repetitive. 
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Careful analysis of the contour plots of the remain- 
ing of the 26 experiments contained essentially the same 


features as those discussed herein. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


The numerical model of Williams (1973) is modified 
in this study to include moisture. The numerical solutions 
become quasi-steady within 1-2 days as in W73. Intensifi- 
cation of the baroclinic zone at mid to upper levels seems 
to be the significant addition to the model due to the in- 
clusion of moisture. This intensification of the baro- 
clinic zone at the mid to upper levels is consistently 
Manifested by the increases of thermal gradient, wind 
shear, and vertical motions in the warm air along the 
frontal zone which are due to latent heat release. Second- 
ly, it should be noted that the model is essentially un- 
changed at low Weavers with the inclusion of moisture. 
Thirdly, the model with moisture shows dependence upon 
the reference potential temperature, G(p,) chosen due 
to the exponential relationship between moisture and tem- 
perature as is indicated by the Clausius-Clapeyron equa- 
tion (3.18). It is recalled that the dry model is 
relatively insensitive to choice of this reference tem- 
perature. The moist model illustrates this exponential 
thermal dependence in the various experiments at the 
BeorometecmOOlLentaat temperatures of 9270, 280, 290, and 
300°K. 

In future studies using Williams model, attention 


should be directed towards transformation to ln p 
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coordinates for a better depiction of upper level gradi- 
ents. Upper level boundaries might also be changed and 
raised to include a stratosphere and finally a more 
sophisticated convective scheme might be applied to the 


model such as that of Arakawa and Schubert (1973). 
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Figure 1. Frontal scale. Experiments 1 and 2. 
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Figure 3. Frontal scale. Experiments 12 and 13. 
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